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ABSTRACT 

This  report  provides  details  of  a  two-component  strain-gauge  balance  and  ancillary  equipment 
that  has  been  developed  to  measure  flow-induced  loads  on  models  in  the  DSTO  water  tunnel. 
The  loads  are  very  small  and  the  balance  was  designed  to  measure  normal  forces  and  pitching 
moments  within  the  ranges  ±2.5  N  and  ±0.02  N.m  respectively.  Due  to  the  small  loads,  it  was 
necessary  to  use  semi-conductor  strain  gauges  on  the  balance.  Balances  having  such  low  load 
ranges  have  not  previously  been  used  within  DSTO.  The  balance  was  used  to  measure  forces 
and  moments  on  a  delta  wing  in  the  water  tunnel  for  a  range  of  fixed  angles  of  attack.  The 
measured  loads  showed  good  agreement  with  wind-  and  water-tunnel  data  reported  in  the 
literature,  showing  that  the  new  load-measurement  system  gives  good  results. 
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Development  of  a  Two-Component  Strain-Gauge-Balanee 
Load-Measurement  System  for  the  DSTO  Water  Tunnel 


Executive  Summary 

Water  tunnels,  sueh  as  the  one  used  at  DSTO,  are  ideally  suited  to  studying  detailed 
flow  patterns  around  models  of  aircraft.  A  shortcoming  of  the  DSTO  water  tunnel,  and 
virtually  all  other  similar  water  tunnels,  is  that  it  was  not  possible  to  measure  forces  and 
moments  on  an  aircraft  corresponding  to  the  observed  flow  patterns.  Simultaneous  load 
measurement  and  visualization  of  the  flow  would  enable  measured  loads  to  be  properly 
interpreted  in  terms  of  the  observed  flow  behaviour.  To  address  this  issue,  a  two- 
component  strain-gauge-balance  load-measurement  system  has  been  developed  for  the 
DSTO  water  tunnel. 

The  main  reason  that  very  few  balances  have  been  developed  for  water  tunnels  is  that 
the  flow-induced  loads  on  an  aircraft  are  small  and  are  difficult  to  measure.  The 
balance  for  the  DSTO  water  tunnel  has  been  designed  to  measure  normal  forces  and 
pitching  moments  within  the  ranges  ±2.5  N  (±0.25  kgf)  and  ±0.02  N.m  (±0.002  kgf  m) 
respectively.  These  loads  are  more  than  1000  times  smaller  than  those  regularly  meas¬ 
ured  on  aircraft  models  in  the  low-speed  wind  tunnel  at  DSTO.  It  was  necessary  to  use 
semi-conductor  strain  gauges  on  the  balance  to  measure  the  small  loads. 

The  balance  was  used  to  measure  normal-force  and  pitching-moment  coefficients  on  a 
delta  wing  in  the  water  tunnel  for  a  range  of  angles  of  attack.  The  measurements  were 
found  to  be  in  good  agreement  with  corresponding  published  wind-  and  water-tunnel 
data,  suggesting  that  the  measured  data  were  credible.  The  capability  of  measuring 
small  flow-induced  loads  on  a  model  in  the  water  tunnel  complements  earlier  work  at 
DSTO  in  which  small  flow-induced  pressures  on  the  surface  of  a  model  in  the  tunnel 
were  measured  for  the  first  time.  The  strain-gauge-balance  load-measurement  system 
has  recently  been  incorporated  into  a  dynamic-testing  system  in  the  water  tunnel, 
enabling  loads  to  be  measured  on  an  aircraft  when  it  is  in  motion.  Details  of  the  strain- 
gauge-balance  load-measurement  system  are  given  in  this  report. 
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Notation 
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[A  ]  Matrix  defined  by  equation  A 1 2 . 

Amt  Ratio  of  delta  wing  projected  cross-sectional  area,  as  viewed  in  the  free- 

stream  direction,  to  test  section  cross-sectional  area,  (non-dimensional). 

C  Coefficient  used  in  the  calibration  equations  for  a  strain-gauge  balance. 

Cm  Pitching-moment  coefficient,  (non-dimensional).  Cm  =  m/(0.5pU  Sc). 

Cn  Normal-force  coefficient,  (non-dimensional),  Cn  =  Z/(0.5pU  S). 

Cw  Empirical  coefficient,  (non-dimensional),  Cw  =  1.57. 

Co  Root  chord  of  delta  wing,  (m),  Co  =  0.3  m. 

c  Mean  aerodynamic  chord  (mac)  of  a  model,  (m). 

\E\  Matrix  defined  by  equation  All. 

e  Sum  of  squares  of  residuals,  (used  when  calculating  calibration  coefficients 

for  a  strain-gauge  balance). 

/  Number  of  degrees  of  freedom  in  the  calibration  equations. 

H  Load  applied  to  a  strain-gauge  balance. 

H  Load  estimated  using  a  calibration  equation. 

L  Linear  correlation  coefficient  (the  correlation  coefficient  indicates  how  well 

a  calculated  curve  fits  the  original  data). 

/,  m,  n  Rolling,  pitching  and  yawing  moments  respectively  associated  with  the 
balance  and  the  model  coordinate  systems,  (N.m).  Positive  directions  are 
given  on  Figures  3  and  9. 

N  Total  number  of  points  used  in  a  calibration. 

P  Index  of  summation. 

R  Electrical  resistance  of  a  strain-gauge,  (Q). 

R  Voltage  ratio,  R  Tout/ Lin* 

S  Plan-view  projected  area  of  a  model  at  a  =  O'*,  (m^). 

se  Standard  error,  used  to  assess  the  accuracy  of  a  calibration. 

U  Free-stream  velocity  in  the  test  section  of  a  tunnel,  (m/s). 

Four,  Fin  Output  and  input  voltages  respectively  for  a  channel  on  a  strain-gauge 
balance,  (V). 

X,  Y,  Z  Axial,  side  and  normal  forces  respectively  associated  with  the  balance  and 
the  model  coordinate  systems,  (N).  Positive  directions  are  given  on 
Figures  3  and  9. 

X,  y,  z  Axes  for  the  balance  and  model  coordinate  systems.  Positive  directions  are 
given  on  Figures  3  and  9. 

Greek  Letters 

a  Angle  of  attack,  (deg). 

ats  Angle  of  attack  corresponding  to  maximum  static  value  of  Cn,  (deg). 

p  Density  of  air  or  water,  (kg/m^). 

Subscripts 

1  Refers  to  the  normal-force  channel  for  a  strain-gauge  balance. 

2  Refers  to  the  pitching-moment  channel  for  a  strain-gauge  balance. 
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1.  Introduction 


The  water  tunnel  at  the  Defenee  Seience  and  Teehnology  Organisation  (DSTO)  has 
been  used  for  many  years  to  study  detailed  flow  patterns  around  models  of  aircraft. 
Water  tunnels  are  better  suited  to  flow-visualization  studies  than  are  wind  tunnels,  due 
to  water  having  a  higher  density  and  lower  mass  diffusivity  than  air,  and  the  fact  that 
free-stream  velocities  used  in  water  tunnels  are  also  substantially  lower  than  those  used 
in  wind  tunnels.  However,  a  shortcoming  of  the  DSTO  water  tunnel,  and  virtually  all 
other  similar  water  tunnels  around  the  world,  is  that  it  has  not  been  possible  to  measure 
forces  and  moments  on  an  aircraft  corresponding  to  the  observed  detailed  flow  patterns 
over  the  aircraft.  Researchers  need  both  the  flow  patterns  and  the  loads  for  different 
aircraft  orientations  to  interpret  how  variations  in  the  flow,  such  as  changing  locations 
for  vortex  breakdown,  lead  to  variations  in  the  loading.  To  address  this  issue,  a  strain- 
gauge-balance  force/moment-measurement  system  has  been  developed  for  the  DSTO 
water  tunnel.  The  features  of  this  new  capability  are  described  in  this  report. 

2.  DSTO  Water  Tunnel 

The  DSTO  water  tunnel  was  manufactured  by  Eidetics  International  Incorporated'  and 
is  designated  Model  1520.  The  tunnel,  shown  in  Figure  1,  holds  about  3790  litres  of 
water  and  has  a  horizontal-flow  test  section  380  mm  wide,  510  mm  deep  and  1630  mm 
long.  It  is  a  recirculating  closed-circuit  tunnel  and  there  is  a  free  water  surface  in  the 
test  section.  The  side  walls  and  floor  of  the  test  section  are  made  from  glass  to  facilitate 
flow-visualization  studies.  The  tunnel  is  constructed  so  that  it  is  possible  to  look 
directly  upstream  into  the  test  section  through  a  glass  window  at  the  downstream  end  of 
the  diffuser.  The  free-stream  velocity  in  the  test  section  can  be  varied  between  0  and 
0.6  m/s.  Downstream  of  the  settling  chamber  there  are  flow-conditioning  elements, 
consisting  of  a  perforated  stainless-steel  plate,  a  honeycomb  and  two  fibreglass  screens. 
The  contraction  upstream  of  the  test  section  has  an  inlet/outlet  area  ratio  of  6:1.  This 
geometry  gives  a  good  velocity  distribution  and  turbulence  level,  and  avoids  the 
likelihood  of  local  flow  separations.  The  free-stream  velocity  in  the  test  section  varies 
by  less  than  ±2%  relative  to  the  average  value,  the  mean  flow  angularity  is  less  than 
±1.0°  in  both  pitch  and  yaw  angle,  and  the  turbulence  intensity  is  less  than  1.0%  RMS, 
(data  as  given  in  the  Eidetics  Water  Tunnel  Operations  Manual). 

There  are  six  dye  canisters  on  the  tunnel  that  can  be  pressurised  with  air  to  force  dye 
through  plastic  tubes  to  selected  locations  on  a  model  for  flow-visualization  studies. 
There  is  a  suction  pump  on  the  tunnel  that  can  be  used  to  suck  or  blow  water  through 
models.  The  suction/blower  system  has  two  independent  circuits,  and  the  maximum 
volume  flow  rates  through  each  circuit  can  be  varied  between  0  and  about  100  ml/s. 
The  flow  rates  through  the  dye  and  suction/blower  circuits  are  controlled  by  hand- 
operated  valves. 

'  Now  called  Rolling  Hills  Research  Corporation,  420  N.  Nash  St.,  El  Segundo,  CA  90245,  USA. 
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Figure  1.  Eidetics  Model  1520  water  tunnel. 

(a)  diagrammatic  view  of  tunnel,  (b)  photograph  of  tunnel. 
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Models  are  mounted  on  a  C-strut  so  that  the  required  centre  of  rotation  of  a  model  is  at 
the  centre  of  the  imaginary  circle  formed  by  the  strut,  ensuring  that  all  angular  motion 
of  the  model  is  about  this  point.  The  model  supporting  system  is  attached  to  the  top  of 
the  test  section  by  a  hinge,  so  that  a  model  can  be  lowered  into  the  test  section,  or 
removed  from  the  test  section,  as  required.  When  a  model  has  been  lowered,  its  centre 
of  rotation  is  at  the  mid  transverse  position  in  the  test  section,  270  mm  from  the  floor 
and  965  mm  downstream  of  the  start  of  the  test  section. 

3.  Strain-Gauge  Balance 

A  sensitive  two-component  strain-gauge  balance,  shown  in  Figure  2,  has  been  designed, 
manufactured  and  tested  to  measure  normal  forces  and  pitching  moments.  The  balance 
is  an  integral  unit,  machined  out  of  stainless  steel  rod,  with  a  specification  AISI  630 
Condition  H480.  The  dimensions  of  the  flexure  member  and  the  likely  deflection  of  the 
balance  under  load  were  determined  using  finite-element  analysis.  A  model  is  bolted 
onto  the  balance  at  its  tapered  end.  The  flow-induced  loads  on  aircraft  in  the  water 
tunnel  are  small  by  conventional  standards,  and  the  balance  has  been  designed  to 
measure  normal  forces  and  pitching  moments  within  the  ranges  ±2.5  N  and  ±0.02  N.m 
respectively.  The  loads  were  determined  by  scaling  force  and  moment  data  measured 
on  full-size  combat  aircraft  during  manoeuvres.  For  comparison,  the  balances  used  in 
the  DSTO  low-speed  wind  tunnel  (LSWT)  are  typically  designed  to  measure  normal 
forces  and  pitching  moments  within  the  ranges  ±3500  N  and  ±300  N.m  respectively. 
The  test  section  of  the  LSWT  has  cross-sectional  dimensions  of  2.74  m  (width)  by 
2.13  m  (height)  and  the  maximum  velocity  obtainable  in  the  test  section  is  about 
100  m/s.  There  is  a  factor  of  1400  for  normal  forces  and  15000  for  pitching  moments 
between  the  loads  for  the  balances  for  the  two  tunnels. 


Figure  2.  Two-component  balance,  showing  the  main  dimensions. 
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C^n 


Gauges: 

1.  Normal  force 

2.  Pitching  moment 


Output 


5V« 


Input 

Normal-force  bridge 


Pitching-moment  bridge 


Figure  3.  Two-component  balance,  showing  the  position  of  the  gauges  for  the  normal- 
force  and  pitching-moment  circuits,  together  with  the  corresponding  bridges. 


Due  to  the  small  loads,  semi-conductor  strain  gauges  have  been  used,  having  a  resist¬ 
ance  of  1000  Q,  a  gauge  factor  of  145,  and  dimensions  of  1.27  mm  by  0.15  mm  (active 
area,  length  by  width).  The  gauges  were  manufactured  by  PSI-TRONIX  Incorporated 
and  the  model  number  of  the  gauges  is  PO 1-05- 1000.  The  balance  contains  four  gauges 
on  each  of  the  two  sides  of  the  flexure  member,  as  shown  in  Figure  3.  Gauges  lb,  Ic, 
2d  and  2b,  which  are  obscured,  are  on  the  opposite  side  of  the  flexure  member  and  are 
directly  underneath  gauges  la.  Id,  2a  and  2c  respectively.  The  eight  gauges  are 
connected  together  to  form  two  Wheatstone  bridges,  as  shown.  The  gauges  have  been 
positioned  on  the  flexure  member  and  wired  together  so  that  the  normal  forces  and 
pitching  moments  measured  by  the  balance  are  not  affected  by  other  load  components, 
such  as  transverse  forces  and  rolling  moments. 


^  PSI-TRONIX  Incorporated,  3950  South  “K”  Street,  Tulare,  CA,  93274,  USA. 
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The  right-handed  orthogonal  balance  coordinate  system  is  shown  in  Figure  3.  The 
coordinate  system  remains  fixed  to  the  balance,  with  the  origin  located  on  the  long¬ 
itudinal  axis  of  the  balance  at  the  geometric  centre  of  the  gauges.  The  directions  of  the 
corresponding  forces  and  moments  are  as  shown,  with  the  positive  directions  of  the 
moments  determined  by  the  right-hand  screw  rule  about  the  respective  axes.  The 
balance  can  only  measure  normal  forces  (Z)  and  pitching  moments  (m). 

The  strain  gauges  have  been  glued  onto  the  balance  using  a  bonding  material  designated 
M-BOND  600  ,  which  is  the  same  as  that  used  for  conventional  strain  gauges.  The 
gauges  and  the  connecting  leads  have  been  waterproofed  by  coating  them  with  a 
compound  designated  M-coat  C^,  a  solvent-thinned  (naptha)  RTV  (room  temperature 
vulcanizing)  silicone  rubber.  Care  was  taken  to  not  use  excessive  amounts  of  the 
waterproofing  agent,  since  this  could  have  stiffened  the  balance  and  caused  it  to  exhibit 
hysteretic  behaviour  during  oscillating  loading.  Figure  4  shows  the  gauges  glued  onto 
the  balance  prior  to  waterproofing.  The  gauges  were  found  to  be  slightly  sensitive  to 
both  temperature  and  light,  but  this  did  not  have  any  significant  effect  on  measured  data. 
By  positioning  the  gauges  back-to-back,  the  effects  of  variations  in  temperature  were 
minimized.  The  temperature  of  the  water  in  the  tunnel  was  very  stable,  varying 
typically  by  less  that  0.2°C  throughout  a  day,  and  the  balance  was  mounted  deep  inside 
a  model  when  testing,  well  protected  from  strong  light  sources  used  when  visualizing 
the  flow. 


Figure  4.  Two-component  balance,  showing  gauges  and  leads,  prior  to  waterproofing. 


^  M-BOND  600  and  M-Coat  C  are  manufactured  by  Vishay. 
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4.  Signal-Conditioning  System 

An  existing  six-channel  signal-conditioning  system  was  modified  and  recommissioned 
for  use  with  the  water-tunnel  balance.  The  system  can  be  used  to  null,  amplify  and  filter 
output  voltages  from  the  Wheatstone  bridges  on  the  balance,  before  the  signals  are 
sampled  by  a  PC-based  data-acquisition  system  (see  Section  5).  A  simplified  circuit 
diagram  of  one  of  the  bridges,  the  power  supplies,  a  nulling  unit,  an  amplifier,  a  filter, 
and  the  data-acquisition  system,  is  shown  diagrammatically  in  Figure  5. 

Bridge  on 

strain-gauge  balance 


C 


Figure  5.  Simplified  circuit  diagram  of  a  bridge  on  the  balance, 
together  with  ancillary  components. 

An  adjustable  power  supply  is  used  to  provide  a  +5.0  V  DC  excitation  voltage  to  the 
strain-gauge  bridges  and  the  nulling  units,  and  another  power  supply  is  used  to  provide 
a  ±15  VDC  excitation  voltage  to  the  amplifiers  and  filters.  Each  channel  contains  a 
bridge  nulling  unit,  with  both  course  and  fine  adjustments,  which  can  be  used  to  null  the 
output  voltage  from  a  Wheatstone  bridge,  i.e.  the  bridge  output  voltage  can  be  set  to 
0  V.  Each  channel  also  contains  an  amplifier,  whose  gain  can  be  set  to  100,  200,  400, 
500  1000  or  2000,  as  well  as  a  third-order  Butterworth  low-pass  filter  whose  3  dB  cut¬ 
off  frequency  can  be  set  to  1,2,  5,  10,  20,  50,  100,  200  or  500  Hz.  The  filter  can  also  be 
switched  out,  in  which  case  the  cut-off  frequency  of  the  filter  is  about  50  kHz,  which  is 
the  inherent  bandwidth  of  an  amplifier. 

When  nulling  a  bridge,  the  voltage  at  the  output  of  a  filter  (points  E  and  F  on  Figure  5) 
is  monitored  rather  than  the  voltage  at  the  output  of  the  bridge  (points  A  and  C  on 
Figure  5).  The  voltage  at  the  output  of  the  filter  can  be  read  from  a  dial  on  a  console,  or 
alternatively,  can  be  read  by  connecting  a  voltmeter  across  the  terminals  labelled 
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OUTPUT  on  the  console.  A  sampling  program,  named  ZERO  BALANCE,  based  on 
LabVIEW"^  software,  has  also  been  developed  so  that  the  voltage  can  be  read  on  a 
monitor  using  a  graphical  user  interface. 

5.  PC  and  Data-Acquisition  Card 

The  PC  used  has  a  Windows™  XP  Professional  operating  system,  a  3.2  GHz  central 
processing  unit,  a  220  GB  hard  disk  drive  and  1  GB  of  RAM. 

The  data-acquisition  card  interfaced  with  the  PC  was  manufactured  by  National 
Instruments™  and  the  product  code  is  NI  6013.  The  card  features  16  channels 
(8  differential)  of  16-bit  analogue  input,  a  68  pin  connector  and  8  lines  of  digital 
input/output.  Analogue  input  voltages  varying  between  -5  V  and  +5  V  can  be  sampled. 
For  the  16-bit  card,  the  resolution  of  the  sampled  voltages  is  10.0/2*^,  i.e.  10.0/65536 
=  0.000153  V/LSB  =  0.153  mV/LSB  (USB  denotes  least  significant  bit). 

When  carrying  out  tests,  it  was  not  possible  to  “see”  the  actual  loads  measured  by  the 
balance.  Loads  had  to  be  determined  from  the  measured  voltages,  using  calibration 
relationships,  after  the  tests  had  been  completed. 

6.  Calibration  of  the  Strain-Gauge  Balance 

6.1  Balance  Calibration  Equations 

The  format  of  the  calibration  relationships  depends  upon  which  particular  balance 
calibration  model  or  set  of  equations  is  chosen  to  represent  the  data.  There  are  a  number 
of  different  sets  of  equations  that  can  be  used,  but  the  one  currently  used  at  DSTO  for 
the  wind-tunnel  balances  assumes  that  balance  output  voltages  are  functions  of  the 
calibration  coefficients  and  the  applied  loads  (see  Lam  1989,  Leung  &  Link  1999, 
Blandford  2004).  This  model  will  also  be  used  for  the  water-tunnel  balance. 

For  a  perfect  balance,  the  output  voltage  from  say  the  normal-force  channel  is  only 
affected  by  the  normal  force  applied  to  the  balance,  and  similarly  for  other  channels. 
However,  in  practice,  balances  are  not  ideal  and  there  may  be  interactions  between  the 
different  channels,  so  that  the  output  voltage  from  each  channel  generally  depends  on  all 
of  the  load  components  applied  to  the  balance. 

For  the  two-component  balance,  the  calibration  is  always  described  by  two  different 
equations,  but  the  number  of  terms  in  each  equation  can  vary  depending  on  whether  the 
equations  are  first-order,  second-order,  third-order  and  so  on.  Both  first-order  and 
second-order  calibrations  were  carried  out  and  the  equations  and  the  calibration 


Lab  VIEW™,  a  graphical  programming  language,  is  a  product  of  National  Instruments. 
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coefficients  are  given  in  this  report,  but  loads  on  a  delta-wing  model  were  computed 
using  first-order  equations  (see  Section  8). 

The  first-order  equations  are  given  by 

(1) 

and 

(2) 

Ri  and  R2  are  voltage  ratios  for  the  normal-force  and  pitching-moment  channels 
respectively,  i.e.  the  output  voltages  from  the  channels  divided  by  the  respective  input 
voltages  to  the  channels,  H\  and  H2  are  the  applied  normal  force  and  pitching  moment 
respectively,  and  the  C  terms  are  calibration  coefficients.  The  second-order  equations 
are  given  by 


(3) 

(4) 

The  two  first-order  equations,  expressed  in  matrix  form,  are  given  by 


and 


R,  =  C,,//,  +  C,,//,  + 


C2  2^2  ^2  \  ^2  12^1  ^2  12^1^2 


'A 

A 

a; 

A. 

A 

Q,2_ 

A2_ 

(5) 


and  similarly,  the  two  second-order  equations,  expressed  in  matrix  form,  are  given  by 


'A 

Cu 

c 

'^1,11 

c 

^\,22 

c 

^^1,12 

A. 

Q,2 

c 

^^2,11 

c 

^2,22 

c 

'^2,12_ 

Equations  5  and  6  can  be  expressed  in  a  simplified  matrix  form  as  follows: 

M=[c][//] 


(6) 


(7) 
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6.2  Acquisition  of  Calibration  Data 

The  balance  was  calibrated  manually  using  the  experimental  setup  shown  in  Figure  6. 
The  balance  was  clamped  onto  a  vee-block  using  a  stirrup  clamp  (not  shown)  so  that  its 
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centerline  was  horizontal  and  so  that  the  positive  direction  of  z  axis  of  the  balance 
coordinate  system  was  vertically  downwards.  The  roll  angle  of  the  balance  was  set  to  0“ 
by  inserting  an  alignment  pin  through  a  hole  in  the  shaft  of  the  balance  and  into  a 
matching  hole  in  the  vee  block.  A  calibration  cage  was  bolted  onto  the  model  end  of  the 
balance  at  the  tapered  adaptor  (see  Figure  3).  The  only  contact  between  the  balance  and 
the  cage  was  at  the  adaptor  and  the  cage  did  not  inhibit  the  bending  of  the  balance.  A 
pan  was  suspended  using  a  knife  edge  along  the  top  of  the  cage  so  that  forces  and 
moments  could  be  applied  to  the  balance  by  placing  weights  in  the  pan.  A  10  mm  by 
10  mm  square  block,  containing  crossed  vees,  was  positioned  between  the  knife  edge  on 
the  cage  and  a  knife  edge  on  the  pan  supporting  arm  to  prevent  damage  to  the  two  knife 
edges.  The  longitudinal  location  of  the  knife  edge  on  the  top  of  the  crossed  vee  block 
determines  the  value  of  x  used  when  computing  applied  moments.  The  maximum 
forces  and  moments  applied  to  the  balance  during  the  calibration  were  about  ±1.4  N  and 
about  ±0.01 1  N.m,  which  covered  the  range  of  loads  likely  to  be  experienced  by  a  delta¬ 
wing  model  during  tests  (see  Section  8). 


Alignment  pin 
Vee  block 


Two-component 
strain-gauge  balance 


Knife  edge  on  cage 

Calibration  cage 
Crossed  vee  block 
level 

Pan 

Weights 


Figure  6.  Experimental  setup  used  when  calibrating  the  two-component  balance. 
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For  the  current  calibration,  the  gain  of  the  amplifiers  was  set  to  500  and  the  cut-off 
frequency  of  the  low-pass  filters  was  set  to  1  Hz.  Prior  to  commencing  each  loading 
sequence,  when  there  were  no  weights  in  the  pan,  the  output  voltages  from  the  two 
channels  were  nulled,  i.e.  the  bridge  output  voltages  were  set  close  to  0.0  V.  The 
method  used  to  null  the  bridges  is  described  in  Section  4. 

Weights  were  progressively  added  to  the  pan,  and  then  progressively  removed,  for  the 
pan  positioned  at  different  x  locations  along  the  knife  edge.  The  x  locations  used  were 
-8.0,  0.0  and  +8.0  mm,  and  for  each  of  these  locations  the  weights  used  were  0  to  140  to 
0  g,  in  increments/decrements  of  20  g.  For  each  of  the  15  loads  used,  the  output 
voltages  from  the  normal-force  channel  and  from  the  pitching-moment  channel  were 
sampled  20  times  at  0.1  s  intervals  and  each  set  of  20  samples  was  then  averaged  to 
obtain  mean  values.  The  balance  was  rolled  180°  and  the  above  loading  sequence  was 
repeated  for  the  same  x  locations,  so  that  the  loading  was  in  the  opposite  direction. 
Altogether,  90  calibration  points  were  taken  for  the  normal-force  channel  and  90 
corresponding  points  for  the  pitching-moment  channel  (15  loads  for  each  of  2  directions 
at  each  of  3  x  locations). 

A  sampling  program,  named  SCAN,  based  on  Lab  VIEW™,  was  developed  to  facilitate 
the  calibration  process.  For  each  of  the  90  pairs  of  calibration  points,  the  user  had  to 
specify  the  number  of  samples,  time  between  samples,  output  file  name,  applied  load, 
excitation  voltage,  temperature,  gain  settings  and  filter  settings.  The  program  was  then 
run  and  the  output  voltages  from  the  different  channels  were  displayed  on  the  monitor  in 
graphical  form  as  they  were  being  sampled.  The  voltages  were  written  to  a  spreadsheet 
file  along  with  other  information  specified  by  the  user.  Input  conditions  were  updated 
for  the  next  pair  of  calibration  points  and  the  procedure  was  repeated.  Six  individual 
data  files  were  created  for  the  complete  calibration. 

6.3  Plots  of  Voltage  Ratios  vs  Loads 

For  the  calibration  data,  plots  of  Four/ Fin  vs  Z  for  the  normal- force  channel  and  plots  of 
Fout/Fin  vs  m  for  the  pitching-moment  channel  are  shown  in  Figure  7.  Z  is  the  applied 
normal  force,  m  is  the  applied  pitching  moment,  Fqut  is  the  output  voltage  from  a 
channel  and  Fm  is  the  input  voltage  to  the  Wheatstone  bridge  for  that  channel  (Fm  is 
nominally  5.0  V).  Each  set  of  calibration  points,  corresponding  to  loading  and 
unloading  the  balance  at  a  given  x  location,  occupies  one  quadrant  on  a  plot.  For 
Figures  7a  and  7c,  the  pitching  moments  are  different  for  each  value  of  Z,  and  for 
Figures  7d  and  7f,  the  normal  forces  are  different  for  each  value  of  m.  Corresponding 
points  for  the  loading  and  the  unloading  phases  are  almost  coincident,  with  no 
discemable  evidence  of  hysteretic  behaviour.  Figure  7e  shows  that  for  pure  normal- 
force  loading  (for  loads  to  ±140  g)  at  x  =  0  mm,  the  output  voltages  from  the  pitching- 
moment  channel  are  very  small,  i.e.  there  is  very  little  cross  coupling,  and  the  30 
calibration  points  are  superimposed.  Where  applicable,  each  set  of  experimental  points 
has  been  curve  fitted  by  a  straight  line  using  the  method  of  least  squares,  and  the 
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equations  for  the  straight  lines  are  as  shown  (R  =  Fqut/I^in)-  The  values  of  the  linear 
correlation  coefficients,  L,  are  also  given.  Each  set  of  experimental  points  closely 
follow  a  straight  line,  as  indicated  by  the  values  of  L  being  close  to  1.0. 

6.4  Evaluation  of  Calibration  Coefficients 

The  calibration  coefficients  were  determined  from  the  discrete  applied  loads  and 
associated  balance  output/input  voltages  using  the  least-squares  regression  method 
proposed  by  Ramaswamy  et  al.  (1987).  Details  of  the  procedure  used  are  given  in 
Appendix  A. 

The  first-order  calibration  coefficients  are 


■Q,! 

Q,2 

“  6.006E-1 

E139E  +  0" 

_Q,1 

Q,2_ 

-4.752E-3 

8.640E  +  1 

(8) 


and  the  second-order  coefficients  are 


■^1.1 

Q,2 

Q,ii 

Q,22 

Q,12 

“  6.006E-1 

E139E  +  0 

E640E-3  5.001E  +  0  ElHE-l" 

_Q,i 

Q,2 

Q,ii 

Q,22 

Q,12_ 

-4.752E-3 

8.640E  +  1 

2.322E-3  -2.087E  +  1  -E070E-1 

6.5  Checking  Accuracy  of  Calibration 

When  calibrating  a  balance,  loads  are  applied  to  the  balance  and  corresponding 
output/input  voltage  ratios  are  measured.  To  check  the  accuracy  of  the  calibration,  the 
reverse  procedure  is  used,  whereby  the  calibration  equations  are  used  to  estimate  loads 
corresponding  to  the  voltage  ratios  measured  during  the  calibration.  For  an  ideal 
calibration,  the  applied  and  estimated  loads  should  be  the  same,  but  in  practice  small 
discrepancies  do  occur.  Normal  forces  and  pitching  moments  were  estimated  for  the 
measured  voltage  ratios  by  using  equation  7.  Since  the  loads  {H  terms)  are  not  an 
explicit  function  of  the  calibration  coefficients  (C  terms)  and  the  voltage  ratios  {R 
terms),  then  the  loads  had  to  be  estimated  iteratively.  Details  of  how  to  do  this  are 
given  by  Fairlie  (1985)  and  Leung  &  Link  (1999).  However,  provided  that  \C\  is  a  non¬ 
singular  square  matrix,  then  [C]'*  (inverse  of  [C])  can  be  formed  and  equation  7  can  be 
rearranged  to  give 


[//]=[c]-'[r]  (10) 

enabling  the  loads  to  be  estimated  directly.  This  equation  can  only  apply  to  the  first- 
order  equations,  and  not  to  the  second-order  equations,  since  [C]  is  a  square  matrix  for 
only  the  first-order  equations. 
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Figure  7a,  b,  c.  Voltage  ratio  vs  Z  for  the  normal-force  channel  for  various  x  locations. 
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m  (N.m) 

Figure  7d,  e,  f  Voltage  ratio  m  for  the  pitching-moment  channel  for  various  x  locations. 
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A  statistical  indicator  commonly  used  to  assess  the  accuracy  of  the  calibration  is  the 
standard  error  (see  Lam  1989  and  Leung  &  Link  1999).  The  parameter  indicates  the 
goodness  of  the  fit  of  the  calibration  equation  and  the  uncertainties  in  the  coefficients 
determined  by  the  least-squares  method.  The  standard  errors  for  the  normal  force  and 
the  pitching  moments  are  given  by  the  following  respective  expressions 


se^  = 


Zkp-«,,p. 


P=1 


N-f 


1 


N-f 


(11) 


(12) 


where  sei  is  the  standard  error  for  the  normal  forces,  //,  p  is  an  applied  normal  force, 
Nj  p  is  the  corresponding  normal  force  estimated  using  the  calibration  equation,  N  is  the 
total  number  of  points  used  in  the  calibration  for  each  channel  (N  =  90),  P  is  an  index  of 
summation  and /is  the  number  of  degrees  of  freedom  in  the  calibration  equations.  (The 
number  of  degrees  of  freedom  is  equal  to  the  number  of  calibration  coefficients  per 
component,  i.e.  /=  2  for  the  first-order  calibration  equations  and  /=  5  for  the  second- 
order  equations).  The  corresponding  terms  for  the  pitching  moments,  i.e.  se2,  and 
p ,  are  similarly  defined.  The  standard  error  given  by  equation  1 1  has  units  of 
Newton  and  that  given  by  equation  12  has  units  of  Newton  meter.  To  convert  the 
standard  errors  to  dimensionless  parameters,  they  are  divided  by  the  appropriate 
maximum  design  load,  either  a  force  or  a  pitching  moment. 

For  the  first-order  calibration,  the  standard  errors  were  found  to  be  0.1%  for  the  normal- 
force  channel  and  0.1%  for  the  pitching-moment  channel.  This  suggests  that  loads 
measured  by  the  water-tunnel  balance  immediately  after  the  calibration  are  of 
acceptable  accuracy. 

6.6  Checking  Stability  of  Calibration 

To  check  the  stability  of  the  calibration  over  time,  the  balance  was  recalibrated  seven 
days  after  the  initial  calibration,  using  the  same  loading  sequence,  and  the  two 
calibrations  were  compared  to  see  if  there  were  any  significant  changes. 

Plots  of  Four/ Fin  VS  Z  for  the  two  calibrations  for  the  normal-force  channel  and  plots  of 
Fout/Fin  vs  m  for  the  two  calibrations  for  the  pitching-moment  channel  are  shown  in 
Figure  8.  Each  set  of  data  points  has  been  curve  fitted  by  a  straight  line  using  the 
method  of  least  squares,  and  the  equations  for  the  straight  lines  are  as  shown. 
Corresponding  data  from  the  two  calibrations  are  almost  coincident,  showing  that  the 
calibration  is  stable  over  the  time  interval  considered. 
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A  reverse  calibration  procedure  was  applied  to  both  calibrations,  whereby  equation  10 
was  used  to  estimate  loads  from  voltage  ratios  measured  during  the  calibrations.  Loads 
were  estimated  using  (a)  the  voltage  ratios  from  the  first  calibration  combined  with  the 
calibration  coefficients  determined  from  the  first  calibration,  and  (b)  the  voltage  ratios 
from  the  first  calibration  combined  with  the  calibration  coefficients  determined  from  the 
second  calibration.  For  the  complete  range  of  calibration  points  considered,  it  was 
found  that  the  variation  in  estimated  normal  forces  for  cases  (a)  and  (b)  was  always 
within  the  range  ±0.8%,  and  the  variation  in  estimated  pitching  moments  for  the  two 
cases  was  always  within  the  range  ±0.3%.  Thus,  there  is  good  agreement  between  the 
first  and  second  calibrations,  confirming  the  above  finding  that  the  calibration  did  not 
change  significantly  over  the  time  interval  considered. 

7.  Checking  Waterproofing  of  Balance 

Tests  were  undertaken  to  see  whether  the  strain  gauges  and  the  electrical  leads  on  the 
balance  were  affected  by  water.  The  balance  was  submerged  in  water  for  a  period  of 
10  hours  between  the  first  and  second  calibrations  (see  Section  6.6),  and  while  under 
water,  the  gauges  on  the  balance  were  active,  i.e.  the  two  bridges  on  the  balance  were 
supplied  with  an  excitation  voltage  of  5  V  DC.  The  two  calibrations  agreed  closely,  as 
shown  in  Section  6.6,  indicating  that  the  waterproofing  agent  on  the  balance  was 
effective. 


8.  Tests  Using  a  Delta  Wing 

The  strain-gauge  balance  was  used  to  measure  forces  and  moments  on  a  70“  delta  wing 
to  see  if  the  results  agreed  with  those  reported  in  the  literature.  A  delta  wing  is  an  ideal 
test  model  to  use  when  comparing  data  obtained  in  different  wind  and  water  tunnels. 
The  reason  for  this  is  that  the  flow  over  a  delta  wing  at  an  angle  of  attack  is  dominated 
by  two  large  bound  vortices  that  are  formed  by  the  rolling  up  of  the  flow  as  it  separates 
along  the  two  leading  edges  of  the  wing.  The  flow  over  such  a  wing  in  different 
facilities  is  consistent  and  well  defined,  and  the  way  the  associated  normal  forces  and 
pitching-moments  vary  with  angle  of  attack  are  well  established. 

8.1  Details  of  Delta  Wing 

The  delta  wing  used  in  the  tests  is  shown  in  Figure  9.  The  wing  is  made  from  clear 
Perspex.  The  leading  edges  of  the  wing  are  square  to  the  leeward  (suction)  surface  for 
the  first  0.5  mm  from  that  surface,  and  they  are  then  bevelled  at  an  angle  of  30°.  For  a 
delta  wing,  the  mean  aerodynamic  chord  (mac)  is  2/3  of  the  root  chord,  i.e.  c  =  (2/3)  Co 
=  200.0  mm.  The  coordinate  system  (x,  y,  z)  is  as  shown.  A  metal  tubular  insert  is 
embedded  in  the  wing  along  its  centreline  to  allow  it  to  be  fitted  onto  a  strain-gauge 
balance,  and  the  insert  is  faired  into  the  windward  and  leeward  sides  of  the  wing  as 
shown.  The  two-component  balance  was  positioned  in  the  insert  so  that  the  origin  of 
the  balance  coordinate  system  was  located  at  the  50%  mac  position  shown  in  Figure  9. 
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Figure  8a,  b,  c.  Voltage  ratio  vs  Z  for  two  calibrations  for  the  normal-force  channel 
for  various  x  locations.  •  first  calibration,  +  second  calibration. 
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m  (N.m) 

Figure  8d,  e,  f  Voltage  ratio  V5  m  for  two  calibrations  for  the  pitching-moment  channel 
for  various  x  locations.  •  first  calibration,  +  second  calibration. 
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Figure  9.  Delta  wing  used  in  the  DSTO  experiments 

8.2  Force/Moment  Tests 

Normal  forces  and  pitching  moments  were  measured  on  the  delta  wing  for  fixed  angles 
of  attack,  a,  varying  from  0°  to  60°,  in  increments  of  2°.  The  data  were  taken  for  a  free- 
stream  velocity,  U,  of  0.1  m/s.  Normal-force  coefficients,  Cn,  are  shown  in  Figure  10, 
and  corresponding  pitching-moment  coefficients.  Cm,  referenced  to  the  30%,  40%  and 
50%  c  positions,  which  are  located  60,  80  and  100  mm  respectively  from  the  leading 
edge  of  the  mean  aerodynamic  chord,  are  shown  in  Figure  11.  The  Cm  data  have  been 
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computed  relative  to  these  three  reference  positions  to  facilitate  comparing  the  data  with 
that  obtained  by  other  researchers.  Cn  and  Cm  were  computed  from  balance  output 
voltages  using  first-order  calibration  equations.  Each  experimental  point  shown  in 
Figures  10  and  11  corresponds  to  50  sets  of  gauge  output  voltages,  sampled  at  1  s 
intervals  per  set. 


Data  were  corrected  for  the  effects  of  tunnel  blockage  using  the  corrections  proposed  by 
Cunningham  &  Bushlow  (1990).  The  corrections  were  only  made  for  values  of  a 
greater  than  those  corresponding  to  maximum  static  values  of  Cn  (as  was  done  by 
Cunningham  &  Bushlow).  The  corrections  were  required  as  a  result  of  a  significant 
expansion  of  the  wing  wake  when  the  wing  stalls.  For  the  current  tests,  this  value  of  a, 
denoted  by  cifts,  was  38°.  Cunningham  &  Bushlow  developed  a  semi-empirical 
relationship,  based  on  the  angle  of  attack  of  a  model  and  its  blockage,  as  follows. 


C  =  C 

N(corrected)  N(uncorrected) 


^w4nt  ^  ) 


(13) 


For  a>  38°,  Cw,  an  empirical  coefficient,  is  equal  to  1.57,  ^mt  is  the  ratio  of  the  model 
projected  cross-sectional  area  (as  viewed  in  the  free-stream  direction)  to  the  test  section 
cross-sectional  area,  and  a '  is  given  by 


f 

a' =  90 

V 


a -a. 


ts 


90  -  a, 


ts  y 


(14) 


Cunningham  &  Bushlow  developed  the  blockage  corrections  for  their  Cn  measure¬ 
ments,  but  they  also  applied  the  corrections  to  their  Cm  measurements  in  a  similar 
manner.  The  corrected  DSTO  Cn  and  Cm  data  are  plotted  in  Figures  10  and  1 1  respec¬ 
tively,  along  with  the  uncorrected  data.  This  method  of  correcting  the  data  was  also 
used  by  Suarez  et  al.  (1994). 

To  help  validate  the  DSTO  Cn  and  Cm  measurements  taken  on  a  stationary  70°  delta 
wing  in  the  water  tunnel,  they  were  compared  with  similar  measurements  taken  on  70° 
delta  wings  in  other  tunnels.  The  corrected  DSTO  Cn  data  are  compared  with  that  of 
Suarez  et  al.  (1994)  (water  tunnel),  Phillis  (1991)  (wind  tunnel),  Cunningham  & 
Bushlow  (1990)  (water  tunnel),  Brandon  &  Shah  (1988)  (wind  tunnel)  and  Wentz 
(1969)  (wind  tunnel)  in  Figure  12.  Corresponding  Cm  data,  computed  for  different 
moment  reference  centres,  is  compared  with  data  for  matching  reference  centres  in 
Figure  13. 

There  is  good  agreement  between  the  DSTO  data  and  that  given  by  other  researchers. 
Variations  between  the  different  sets  of  data  could  be  due  to  differences  in  Reynolds 
numbers,  flow  quality,  flow  blockage,  and  the  fact  that  the  leading  edges  of  the  wings 
were  not  always  the  same.  For  example,  the  leading  edges  of  the  DSTO  delta  wing 
were  only  bevelled  on  one  side  (see  Figure  9),  whereas  the  leading  edges  of  the  delta 
wing  used  by  Suarez  et  al.  (1994)  were  bevelled  on  both  sides. 
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Figure  10.  Normal-force  coefficients  for  the  DSTO  delta  wing. 
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Figure  11.  Pitching-moment  coefficients  for  the  DSTO  delta  wing. 
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Figure  12.  Normal-force  coefficients  for  the  DSTO  delta  wing  compared  with  other  data. 


Angle  of  Attack  a  (deg) 


Figure  13.  Pitching-moment  coefficients  for  the  DSTO  delta  wing  compared  with  other  data. 


21 


DSTO-TR-1835 


9.  Concluding  Remarks 

This  report  provides  details  of  a  strain-gauge-balance  force/moment  measurement 
system  that  was  developed  to  measure  flow-induced  loads  on  models  in  the  DSTO 
water  tunnel.  The  loads  are  small  by  conventional  standards  and  the  balance  was 
designed  to  measure  normal  forces  and  pitching  moments  within  the  ranges  ±2.5  N  and 
±0.02  N.m  respectively.  Strain-gauge  balances  having  such  low  load  ranges  have  not 
been  used  previously  within  DSTO.  For  comparison,  the  balances  used  in  the  low- 
speed  wind  tunnel  at  DSTO  regularly  measure  loads  greater  than  1000  times  these 
values.  To  measure  the  small  loads,  it  was  necessary  to  use  semi-conductor  strain 
gauges  on  the  balance.  The  calibration  of  the  balance  was  found  to  be  accurate  to 
within  ±0.1%  and  also  stable,  and  the  waterproofing  of  the  balance  was  found  to  be 
effective. 

The  new  system  was  used  to  measure  loads  on  a  70°  delta  wing  for  a  range  of  static 
angles  of  attack.  Measured  normal-force  and  pitching-moment  coefficients  agreed  well 
with  wind-  and  water-tunnel  data  reported  in  the  literature,  showing  that  the  new  system 
gives  good  results,  at  least  for  a  sharp-edged  delta  wing,  where  flow  patterns  are 
independent  of  Reynolds  number. 

The  ability  to  measure  the  small  flow-induced  loads  on  a  model  in  the  water  tunnel 
increases  the  usefulness  of  the  tunnel,  and  complements  earlier  work  done  at  DSTO,  in 
which  the  small  flow-induced  pressures  on  the  surface  of  a  model  in  the  tunnel  were 
measured  (Erm  2000).  The  load-measurement  system  has  recently  been  incorporated 
into  a  dynamic-testing  system  that  enables  loads  to  be  measured  on  an  aircraft  model  as 
it  undergoes  a  defined  dynamic  manoeuvre  (Erm  2006). 
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Appendix  A:  Evaluation  of  Calibration  Coefficients 

In  Section  6.1,  it  was  indicated  that  the  first-order  calibration  relationships  for  the  two- 
component  balance  were 


^1  -  ^1,1^1  +  C'l,2^2 

and 

/?2  ~  +  C2  2^2 


and  the  corresponding  second-order  calibration  relationships  were 


and 


R,  =  +  C,22m  + 


~  C2  1^1  ^2  1^1  ^2  ^2  22^2  ^2  12^1^2 


(Al) 

(A2) 


(A3) 

(A4) 


R\  and  R2  are  voltage  ratios,  viz.  the  output  voltages  from  the  normal- force  and 
pitching-moment  channels  divided  by  the  input  voltages  to  those  channels  respectively, 
H\  and  Hi  are  the  applied  normal  force  and  pitching  moment  respectively,  and  the  C 
terms  are  calibration  coefficients. 


The  coefficients  appearing  in  these  equations  can  be  determined  from  the  discrete 
applied  loads  and  associated  voltage  ratios  using  a  least-squares  regression  model. 
Various  types  of  regression  methods  have  been  reported  in  the  literature,  but  the  method 
used  when  calibrating  the  wind-tunnel  balances  in  Air  Vehicles  Division  is  that 
developed  by  Ramaswamy  et  al.  (1987).  This  method  was  also  used  when  calibrating 
the  water-tunnel  balance.  Lam  (1989)  gives  mathematical  details  of  the  method,  and 
only  an  outline  of  the  method  will  be  given  here  -also  see  Leung  &  Link  (1999)  and 
Blandford  (2004).  The  method  will  only  be  described  for  the  second-order  calibration 
equations  (equations  A3  and  A4),  but  the  same  principles  apply  for  equations  of  lower 
or  higher  orders. 

For  this  method,  the  calibration  coefficients  are  established  when  the  sum  of  the  squares 
of  the  differences  between  the  measured  strain-gauge  output  voltage  ratios  and  those 
obtained  from  a  calibration  equation  are  a  minimum.  That  is. 


and 


^1  -  +  Q,2-^2,P  +  +  ^1.22^1?  +  Q, 12-^1, p-^2, P  “-^l.p] 


= el.  jifj  p  +  C2  2H2  P  -|-C2  j[/7j^p  C2,22^2,P  '^^2,n^\,?^2,?  ~ -^,P  ] 


(A5) 

(A6) 
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are  both  a  minimum,  where  e\  and  ej,  are  sums  of  squares  of  residuals  for  the  normal- 
force  and  pitching-moment  components  respectively,  P  is  an  index  of  summation  and  N 
is  the  number  of  calibration  points  for  each  channel  (N  =  90  -see  Section  6.2).  The 
procedure  involves  partial  differentiating  e\  with  respect  to  each  of  its  coefficients  and 
equating  the  resultant  expressions  to  zero,  and  likewise  for  e2-  This  results  in  a  set  of  5 
non-linear  simultaneous  equations  for  each  of  the  two  components  as  follows^: 


and 


y  +  C,  ,,//,■  +  +  C,  ]//,=  0 

=0 

,=  +  - i!,] ff,’  =  0  (A7) 

ZkuW,  +  C,,H,  +  ,=  +  - 7!,] =  0 

+  C,,,«  +  -  A,]  H,H,_  =  0 


'Z[C,,H,  +  C,,H,  +  -R,]h,=0 

y  [C„//,  +  C„H,_  +  C,  „H;  +  Q  -R,]h,=0 

+  +  Q -R^]hI=0  (A8) 

y  [C„//,  +  C„//,  +  Q  ,,//■  +  C,  „Hl  +  C  „//,//,  -  //’  =  0 

X[C,,7/,  +  +  C  +  C„Ml  +  C, - R-]h,H,  =  0 


Equations  A7  and  A8  can  be  combined  and  put  into  matrix  notation  as  follows: 


[£][c]^=M 


[C]^  is  the  transpose  of  the  matrix  of  calibration  coefficients  and  is  given  by 


(A9) 


^  For  clarity,  the  subscript  p  on  i?  and  H  is  omitted  in  the  remainder  of  this  report.  Whenever  the 
summation  symbol  is  used,  it  is  to  understood  that  the  sum  is  to  be  taken  over  the  range  P  =  1  to  N. 


26 


DSTO-TR-1835 


c  c 

'^l,!  ^^2,1 

[cr= 

c  c 

^\,2  ^^2,2 

c  c 

'^1,11  '^2,11 

c  c 

'^1,22  '^2,22 

c  c 

_'^1,12  '^2,12  _ 

[E]  and  [A]  are 

given  by 

^H,Hl 

[£]= 

T.HIHI 

T.HIHI 

ILHlHl 

(AlO) 


(All) 


and 


Uh 


TH^R^ 

Y.HlR^ 

ILHlR^ 

(A12) 


The  25  elements  eomprising  [E]  and  the  10  elements  comprising  [A]  can  each  be 
determined  from  the  known  applied  calibration  loads  and  the  corresponding  measured 
voltage  ratios,  so  that  the  10  elements  in  [C]^  can  be  determined  by  making  use  of 
equation  A9.  It  is  not  possible  to  determine  [C]^  directly  from  this  equation,  and  the 
equation  must  be  rearranged  to  obtain  an  explicit  expression  for  [C]^.  Provided  [E]  is  a 
non-singular  square  matrix,  then  [£]’'  (inverse  of  [£■])  can  be  formed  and  equation  A9 
can  be  rearranged  as  follows: 


[c]’=[£]-'p]  (A13) 

Once  [C]^  is  known,  [C]  can  be  determined  by  simply  transposing  [C]^  (equation  AlO) 
to  give 


[C] 


c  c  c  c 

'^l,!  ^^1,2  '^l,!!  ^^1,22 

c  c  c  c 

^2,1  ^2,2  ^2,11  ^2,22 


Q,. 

C. 


2,12 


(A14) 
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